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(Dated: August 2, 2008) 

Wc describe a search for Z boson pair production in pp collisions at \/s= 1.96 TcV with the 
DO detector at the Fermilab Tevatron Collider using a data sample corresponding to an integrated 
luminosity of 2.7 fb"^. Using the final state decay ZZ — » t^t~vv (where £ = e or /x) we find a 
signal with a 2.6 standard deviations significance (2.0 expected) corresponding to a cross section of 
cr(pp ZZ + X) = 2.01 ± 0.93(stat.) ± 0.29(sys.) pb. 



PACS numbers; 12.15.-y, 12.15.Ji, 14.70.Hp, 07.05.Kf, 29.85.Fj. 
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INTRODUCTION 

We report a search for Z boson pair production in 
collisions in the mode where one Z boson decays into 
two charged leptons (either electrons or muons) and the 
other Z boson decays into two neutrinos (see Fig.[T]). In 
the Standard Model (SM), ZZ production is the double 
gauge boson process with the lowest cross section and 
is the last remaining unobserved diboson process at the 
Tevatron, aside from the expected associated production 
of the Higgs boson. Observation of ZZ production there- 
fore represents an essential step in Higgs boson searches 
in the ZH and WH channels with sensitivity at the level 
of the expected SM cross sections. Additionally, the ZZ 
process forms a background to Higgs boson searches, for 
example in the channels ZH i^i~bb, ZH i/Vbb and 
H W+W- l+v(.-V. Unhke the WW and WZ 
processes, there are no expected SM contributions from 
triple gauge boson couplings involving two Z bosons and 
a measurement of the ZZ cross section represents a test 
for production of this final state via anomalous couplings. 

The process ZZ I'^t^vv has a branching ratio six 
times larger than that for the other purely leptonic pro- 
cess ZZ l~ l'^ (!~ . After removing instrumental 
backgrounds, the dominant background in the ZZ 
t^l~vv search arises from the process WW V^vl^v, 
which produces the same final state particles. A kine- 
matic discriminant against background from WW — > 
l'^vl~v is employed. In contrast, a search in the ZZ 
t^l~ channel benefits from having no significant 
backgrounds from physics processes with the same final 
state. 

We select events containing an electron or muon pair 
with high invariant mass and significant missing trans- 
verse momentum. After the initial event selection, the 
dominant source of instrumental background to this sig- 
nature arises from events containing a leptonic Z boson 
decay in which the apparent missing transverse momen- 
tum arises from mismeasurement of the transverse mo- 
mentum of either the charged leptons or the hadronic 
recoil system. We introduce a variable that is highly dis- 
criminating against such instrumental background. 

Although the e+e^ ZZ process has been observed 
at LEP p|, ZZ production has not yet been observed 
at a hadron collider where different physics processes are 
allowed and higher energies can be probed. The DO col- 
laboration has previously performed a search for the pro- 
cess ZZ £+£-£'+£'- with £' = e or 0, which set a 
limit on the cross section of a{ZZ) < 4.4 pb and also ex- 
amined non SM triple gauge boson couplings. The CDF 
collaboration has recently produced a result using both 
the £+£-£'+£'- and the £+£-^iyD channels measuring 
the cross section to be a{ZZ) — 1.4^Q g pb. 

The DO detector [3, IE S| contains tracking, calorime- 
ter and muon subdetector systems. Silicon microstrip 




FIG. 1: Leading order Feynman diagram for the process 
ZZ ^e+£-uiy. 



tracking detectors (SMT) near the interaction point cover 
pseudorapidity I77I < 3 to provide tracking and vertexing 
information. The SMT contains cylindrical barrel layers 
aligned with their axes parallel to the beams and disk 
segments. The disks are perpendicular to the beam axis, 
interleaved with, and extending beyond, the barrels. The 
central fiber tracker (CFT) surrounds the SMT, provid- 
ing coverage to about {\r]\ = 2). The CFT has eight 
concentric cylindrical layers of overlapped scintillating 
fibers providing axial and stereo (±3°) measurements. A 
2 T solenoid surrounds these tracking detectors. Three 
uranium-liquid argon calorimeters measure particle en- 
ergies. The central calorimeter (CC) covers jryj < 1, 
and two end calorimeters (EC) extend coverage to about 
I77I — 4. The calorimeter is highly segmented along the 
particle direction, with four electromagnetic (EM) and 
four to five hadronic sections, and transverse to the par- 
ticle direction with typically A77 = A(j) = 0.1, where is 
the azimuthal angle. The calorimeters are supplemented 
with central and forward scintillating strip preshower de- 
tectors (CPS and EPS) located in front of the CC and 
EC. Intercryostat detectors (ICD) provide added sam- 
pling in the region 1.1 < |?7| < 1.4 where the CC and EC 
cryostat walls degrade the calorimeter energy resolution. 
Muons are measured with stations which use scintillation 
counters and several layers of tracking chambers over the 
range \ri\ < 2. One such station is located just outside the 
calorimeters, with two more outside 1.8 T iron toroidal 
magnets. Scintillators surrounding the exiting beams al- 
low determination of the luminosity. A three level trigger 
system selects events for data logging at about 100 Hz. 
The first level trigger (LI) is based on fast custom logic 
for several subdetectors and is capable of making deci- 
sions for each beam crossing. The second level trigger 
(L2) makes microprocessor based decisions using multi- 
detector information. The third level trigger (L3) uses 
fully digitized outputs from all detectors to refine the 
decision and select events for offline processing. 
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DATA SET AND INITIAL EVENT SELECTION 

The data for this analysis were collected with the DO 
detector at the Fermilab Tevatron pp Collider at a center- 
of-mass energy ^/s = 1.96 TeV. An integrated luminosity 
of 2.7 fb^^ is used after applying data quality require- 
ments. The data are selected using a combination of 
single electron or single muon triggers for the respective 
dilepton channels. 

The data taking period prior to March 2006 is refered 
to as Run Ila, while lib denotes the period after. This 
division corresponds to the installation of an additional 
silicon vertex detector, trigger upgrades, and a significant 
increase in the rate of delivered luminosity. 

In each of the two channels we require that there be 
exactly two oppositely charged leptons with transverse 
momentum pT > 15 GeV and dilepton invariant mass 
70 < Mu < 110 GeV. Electrons are required to be within 
the central (I77I < 1.1) or forward (1.5 < I77I < 2.5) regions 
of the calorimeter. We do not use electron candidates 
which point towards the transition region of the central 
and forward cryostats. Electrons must pass tight selec- 
tion criteria on the energy close to them in AR, where 
AR is the distance between two objects in (77, 0) space, 
AR — ^ {^(p)^ + (i5^)^, by requiring that 

^tot(0.4)-gEM(0.2) 

i?EM(0.2) < 

where £"101(0.4) is the total energy within a cone of 
AR < 0.4 and iJ^M (0.2) is the EM energy within a cone 
of AR < 0.2. Additionally, a seven parameter multi- 
variate discriminator compares the energy deposited in 
each layer of the calorimeter and the total shower energy 
to distributions determined from electron GEANT Monte 
Carlo (MC) simulations This discriminator also uses 
the correlations between the various energy distributions 
to ensure that the shower shape is consistent with that 
produced by an electron. 

Each muon is required to have an associated track in 
the central tracking system which has at least one hit 
in the SMT and a distance of closest approach to the 
primary vertex in the plane transverse to the beam of 
|5| < 0.02 cm. Furthermore, the muons must be isolated 
in both the calorimeter and the tracker. For the former, a 
requirement is made that the sum of calorimeter energies 
in cells within an annulus 0.1 < AR < 0.4 around the 
muon track is smaller than 10% of the muon px'- 

E™"^£t(0.1 <AR< 0.4)/pt(m) < 0.1. 

For the latter, the sum of track pT within a cone AR < 
0.5 around the muon track (not included in the sum) 
must be smaller than 10% of the muon px- 

E*™"V(Ai? < 0.5)/pT(/i) < 0.1. 



To suppress background from WZ production, we veto 
events with one or more leptons (e, /Lt, or r) in addition 
to those forming the Z candidate. Additional lepton can- 
didates must be separated by AR > 0.2 from both Z- 
candidate leptons. Electron candidates used in the veto 
must have Et > 5 GeV and either a central track match 
or satisfy shower shape requirements. Muons are rejected 
based on looser quality requirements than those from Z 
decay. Multi-prong hadronic taus are used to form the 
veto if they have been identified using the standard DO 
algorithms [8|. Finally, events are vetoed if they have 
any isolated tracks with pT > 5 GeV and a separation 
distance between the track intercept with the beam line 
and the primary vertex satisfying |Az| < 1 cm. 

Events with relatively large calorimeter activity are re- 
jected by vetoing on the presence of more than two jets 
in the detector. These jets are reconstructed using the 
Runn Ila cone algorithm Q with a radius of 0.5 and 
must satisfy Ai?(jet, lepton)> 0.3 and jet Et > 15 GeV. 
This requirement significantly reduces background from 
tt production. 

Missing transverse energy (^t) is the magnitude of the 
vector sum of transverse energy above a set threshold in 
the calorimeter cells, corrected for the jets in the event. 
At the dilepton selection stage, we do not make a require- 
ment on ■ 

BACKGROUND AND SIGNAL PREDICTION 

Background yields were estimated using a combina- 
tion of control data samples and MC simulation. The 
primary background after the initial selection is inclu- 
sive Z/j* — ^ production. After making the final 
selection described later, the dominant background is 
W~^W^ t'^vd^v events. Additional backgrounds in- 
clude ti production, WZ production and or I^-l-jets 
events in which the 7 or jet is misidentified as an electron. 

The WW, WZ, Z/7*, and tt backgrounds are esti- 
mated based on simulations using the pythia [ld| event 
generator, with the leading order CTEQ6L1 [ll| pa- 
rameterization used for the parton distribution functions 
(PDFs). We pass the simulated events through a de- 
tailed DO detector simulation based on GEANT and re- 
construct them using the same software program used to 
reconstruct the collider data. The Z/7* MC events are 
assigned a weight as a function of generator level jor, to 
match the px spectrum observed in unfolded data [12| . 
Randomly triggered collider data events are added to the 
simulated pythia events. These data events are taken 
at various instantaneous luminosities to provide a more 
accurate modeling of effects related to the presence of 
additional pp interactions and detector noise. We also 
apply corrections for trigger efficiency, reconstruction ef- 
ficiency and identification efficiency. The corrections are 
derived from comparisons of control data samples with 
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simulation. 

The 14^7 background is estimated from a calculation of 
the Next to Leading Order (NLO) production cross sec- 
tion [isj , which we use to normalize events generated by 
PYTHIA. The probability for a photon to be misidentified 
as an electron is measured in a Z — > £^£^j control data 
sample. This probability is then applied to the W'^ yield 
predicted from simulation to determine the contribution 
to our selected sample in which the photon is mistakenly 
reconstructed as the second electron. 

The kinematic distributions of the VF+jet events are 
determined from simulation based on ALPGEN IJ]. The 
overall yield from M^+jet events is determined from data. 
The probabilities that an electron or a jet which satisfy 
looser selection requirements will also pass our candidate 
selection in data are measured by solving a set of lin- 
ear equations involving these probabilities, the number 
of candidate events, and the number of events in which 
one of the requirements on one of the candidates has been 
loosened. The solution to this set of linear equations is 
used to determine the number of W-t-jet events in the fi- 
nal sample. The fraction of jets misidentified as electrons 
is 2 X 10"'^, and the probability that a jet fakes a muon 
is 4 X 10^'' (averaging over Run Ila and Run lib). 

The relative normalization of the background sources 
determined from simulation is taken from ratios of NLO 
cross sections, and the absolute normalization of the total 
background is then determined by matching the observed 
yield under the Z dilepton mass peak to the predicted 
background sum after applying the basic Z selection and 
extra-activity event veto. We choose this normalization 
method because inclusive Z production dominates our 
signal by four orders of magnitude, and this approach 
allows cancellation of multiplicative scale uncertainties 
and systematics related to the modeling of the selection 
efficiencies. The normalization factor agrees with that 
obtained using the integrated luminosity to within the 
associated 6% uncertainty. 

Signal events were generated using pythia with the 
CTEQ6L1 PDFs, and the signal event samples were cor- 
rected for the same detector effects as the background 
samples. 



than 10^, stringent selection criteria against inclusive Z 
production are required. We present a novel approach to 
this challenge. In particular, we do not attempt to make 
an unbiased or accurate estimate of the missing trans- 
verse momentum in the candidate events. Rather, our 
approach is to construct a variable which is a rep- 
resentation of the minimum feasible given the mea- 
surement uncertainties of the leptons and the hadronic 
recoil. Thus, this is not intended to be the best estimator 
of true ^T, but rather to be robust against reconstruc- 
tion mistakes. This approach is inspired by the OPAL 
collaboration which used a similar variable to search in 
final states similar to that of this analysis [HI . 
The variable is constructed in five steps. 

i. The first step is the computation of a reference axis 
chosen such that effects from leptonic resolution oc- 
cur dominantly along this axis, and the decomposi- 
tion of the dilepton system transverse momentum 
into components parallel and perpendicular to this 
reference axis. 

ii. The second step is determination of a recoil vari- 
able based on measured calorimeter jet or total 
calorimeter activity. In events with no significant 
neutrino energy or mismeasurement, the quantities 
calculated in the first two steps should be approxi- 
mately balanced. 

iii. The third step is the calculation of a correction 
based on recoil track px for tracks which are well 
separated from the candidate leptons and calorime- 
ter jets. 

iv. The fourth step is the computation of a correction 
term accounting for lepton transverse momentum 
measurement uncertainties. 

V. The final step is a combination of the quantities 
computed in the first four steps into the vari- 
able. 



Decomposition 



VERIFICATION OF MISSING TRANSVERSE 
MOMENTUM 

The basic event signature for this analysis is a high 
mass pair of charged leptons from the decay of a Z boson, 
produced in association with significant missing trans- 
verse momentum arising from the neutrinos produced 
in the decay of a second Z boson. Substantial background 
comes from single inclusive Z production in which mis- 
measurement results in a mistakenly large value. Be- 
cause the cross section times branching ratios for inclu- 
sive Z production and ZZ signal differ by a factor of more 



In the first step, to minimize the sensitivity to mis- 
measurement of the pt of the individual leptons, the px 
of the lepton pair is decomposed into two components, 
one of which is almost insensitive to lepton px resolu- 
tion for Z candidates with moderate values of transverse 
momentum. This decomposition is achieved as follows. 
In the transverse plane a dilepton thrust axis is defined 
(see Fig. [2]). This axis maximizes the scalar sum of the 
projections of the px of the two leptons onto the axis. It 
is defined as 

r_ ^ (1) _ ^ (2) 

I — Pj, Prp , 
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thrust axis 




FIG. 2: Representation of the transverse plane of the event 
and of the decomposition of the dilepton transverse momen- 
tum along the thrust axis. 

where prj^^^ and p^^^ are the transverse momenta of the 
higher and lower pT leptons respectively. 

We then define two unit vectors di and dt which are 
parallel and perpendicular to the thrust axis. For the rest 
of this paper, a subscript I denotes the component in the 
di directio and a subscript t denotes the component of a 
vector in the dt direction. 

The dilepton system transverse momentum is decom- 
posed into components parallel to di (aj^) and perpen- 
dicular to di (a"). These are given by 

a" = Pt ■ at 
a" = Pt -ai, 

in which pg?^ = p^^' Pt^^ is the dilepton system trans- 
verse momentum. The resolutions of the two components 
are shown in Fig.[3]from Z — + ///i (MC) generated events. 
Resolution effects are more pronounced in the Z — > /i/i 
channel than in the Z — > ee. As seen, the lepton momen- 
tum resolution effects are significantly more pronounced 
in the di direction than in the dt direction. 

The decomposition is performed only for events in 
which A0^^ > 7r/2, where A^^^ is the angle between the 
two charged leptons in the transverse plane. For the case 
A^^^ < 7r/2 the direction of the dilepton transverse mo- 
mentum ^ is used to define dt , and all components in 
the di direction are set to zero. 

Calorimeter Recoil Activity 

The second step in the process uses calorimeter en- 
ergy to assess whether or not it might generate appar- 




a, difference (GeV) a, difference (GeV) 



FIG. 3: The difference between measured and true dilepton 
Pt in simulated Z — > /i/i events projected along the (a) di and 
(b) dt axes. 

ent, though false, I^t- Two measures of net calorimeter 
transverse energy are considered: (a) a vector sum of 
the Et of selected reconstructed jets and (b) the uncor- 
rected missing transverse energy ^t for the event from 
the calorimeter. When computing the jet pt sum, we 
consider only those jets whose pt is in the direction op- 
posite to the dilepton system for each of the di and dt 
directions. That is, if p^^*^* • a; < 0, then this amount is 
added to the di correction, and if pt^*^* - dt < 0, it is added 
to the dt correction. 

The calorimeter activity correction is then defined us- 
ing either the jets or the uncorrected I^t- We chose the 
one with the largest projected magnitude along each of 
the two axes in the hemisphere opposite to the dilepton 
pair. Additionally, we allow for the possibility that only 
some of the true recoiling energy is underestimated by 
multiplying the observed energy by two. The calorime- 
ter correction is thus defined as 

fo™' = 2 X min(EpT^^*' • a*, -fr ■ dt, 0) 
(5a^°' = 2 X min(I]pT^^*' • a;, -^t • a;, 0) 

in which a jet is used in the di{dt) component sum only 
if it satisfies pt ■ di < {pt ■ dt < 0). The presence 
of the zero term in the min-function ensures that the 
calorimeter activity correction is used only if it decreases 
the apparent value of af^ and/or af^. In this way we try 
to minimize the possibility that a well-balanced event ac- 
quires an apparently significant net missing momentum 
due to calorimeter noise or a jet with a grossly overesti- 
mated energy. 

Recoiling Tracks 

The third step identifies events in which the recoil ac- 
tivity is not observed in the calorimeter as jets. We con- 
sider tracks that are AR > 0.5 away from all calorimeter 
jets, AT? > 0.5 from the candidate leptons, have a fit sat- 
isfying x^/NDF < 4.0, and pt > 0.5 GeV and use these 
to build track jets using a cone algorithm. A track jet is 
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seeded by the highest px track not yet associated with 
any track jet. AU tracks which are within AR < 0.5 of the 
seed track and are not yet associated with a track jet are 
added to the current track jet. The track jet transverse 
momentum p^^* is the vector sum of the px values of all 
tracks forming the track jet. This is repeated until no 
unused tracks outside the calorimeter jet cones remain. 
A track-based correction is then defined as 

As with the calorimeter jets, a track jet is included in 
the correction for the di{dt) direction only if it satisfies 
p^^"'^ ■ di < [px"*' • a* < 0). 



Lepton pt Uncertainty 

In the fourth step of the algorithm, corrections (5af^ and 
(5af^ arising from the uncertainties in the lepton trans- 
verse momenta are derived. The basic approach taken is 
to fluctuate the lepton transverse momenta by one stan- 
dard deviation of their uncertainty so as to minimize, 
separately, the di and dt components of the dilepton p^ . 
The transverse component, af^, is minimized by decreas- 
ing the transverse momenta of both leptons to give the 
modified quantity: 

«' ^ «' -I 
^Pt ■ at ■ 

Here pr^^ and dt correspond, respectively, to p^^ and 

dt, redefined using p^j;^"' x (1 ^ ci) and prj?^ x (1 ^ (^2) in 
place of the unsealed quantities. The uncertainty is then 
simply given by: 

The longitudinal component, a^^, is minimized by de- 
creasing p^^^ and increasing p,^^^ using their fractional 
uncertainties cti and (T2- 

If the fractional uncertainty on either of the lepton trans- 
verse momenta is larger than unity, then the fractional 
uncertainties on both af^ and af^ are set to unity. 

Electrons falling at calorimeter module boundaries of 
the central calorimeter require special treatment as their 
calculated uncertainties do not reflect the probability for 
such electrons to have very significantly underestimated 
energies. To account for this, if the lower pr electron is 
within a central calorimeter module boundary, then the 
fractional uncertainty on is set to unity. 



Combination 

In the fifth and final step, the variable is computed 
from the quantities calculated in the previous steps. We 
compute components: 

at = af + (5ar' + fc' x ^^T^ + kx ba\^ 
ai = af -f^ Ja^"' + fc' x Saf'' + kx (5af , 

where k and k' are constants defined below. Recall that 
by construction the Sai (where i = t or I) terms are 
always zero or negative while af^ is positive. 

For events with significant transverse energy from neu- 
trinos and no mismeasurements, the variables are large 
and positive. If ai < 0, then there is no significant miss- 
ing transverse momentum along direction i and that com- 
ponent is ignored in the subsequent analysis by setting: 

a't = max(af,0), 
a'l — max(a;,0). 

The final discriminating variable is then calculated as 
a weighted quadrature sum of the two components: 

^T = \/«f + (l-5a;)2. 

By construction is less that ]f,T- The factor of 1.5 
is used with the dt component to give extra weight to 
the better-measured direction. As mentioned earlier, if 
A</)^^ < 7r/2, then instead of using the thrust axis, the 
reference axis direction dt is simply the p-/^ direction, 
and the a; components are ignored. 

The values k and fc' were optimized by applying a loose 
cut in J/^'rj. (such that the background in the sample is 
dominated by Z events) and maximizing Sj ypB, 

where S is the number of signal events and B is the num- 
ber of background events. The chosen values for dielec- 
tron events are fc = 2.2 and fc' = 2.5. For dimuon events, 
k = k' = 1.5. 

The power of the variable is displayed in Figure IH 
which shows the distribution of for the dielectron 
and dimuon channels separately. 

The separation of sources with true (e.g. WW and 
ZZ) from those without, especially inclusive Z produc- 
tion, is clearly visible. The rejection of events from single 
Z boson production using this method and a simple ]/^t 
cut is shown in Figure [5l 

FINAL SELECTION, LIKELIHOOD AND YIELDS 

In addition to the initial selection requirements, events 
selected for further analysis must satisfy, 

^'j, > 27 GeV, dielectron, Ila, 

> 30 GeV, dimuon, Ila, 

> 27 GeV, dielectron, lib, 

> 35 GeV, dimuon, lib. 
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FIG. 4; The distributions for dielectron (a) and dimuon 
(b) events. 
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FIG. 5: The efficiency of WW, WZ, and ZZ events vs single 
Z events using the (labeled as MET in the figure) and 
variables in the dielectron (a) and dimuon (b) channels. 



The values of these requirements were chosen so as to 
optimize the expected significance of a ZZ observation 
in the four individual analysis channels, assuming the 
SM cross sections for signal and backgrounds. The ef- 
fect of systematic uncertainties, as described below, were 



used in the significance calculation for the optimization. 
Tables U and [IT] show the predicted and observed yields 
after the initial selection and after the ^'rp selection for 
the dielectron and dimuon channels respectively. The 
requirement on ^'j, reduces the predicted instrumental 
background yields well below those for our signal and re- 
maining physics backgrounds. In the dielectron channel, 
we observe 8 events (8.9 expected) in the Ila data, with 
another 20 events (10.7 expected) in the lib data. Of 
these, we expect 1.8 and 2.3 to be signal events respec- 
tively. In the dimuon channel, we observe 10 events (7.0 
expected) in the Ila data and 5 events (7.3 expected) in 
the lib data. Here, we expect 1.7 signal events in each 
data set. 



TABLE I: Number of predicted dielectron events and yield 
observed in data after the dilepton selection and after the 
requirement on 1?^ . The uncertainties in the final column are 
statistical only. If not present, the statistical uncertainty is 
negligible. 



Sample 


dilepton selection 


]^'rp requirement 


Z £+£- 


1.18 X 10'' 


0.5 ± 0.2 


Z -* T+T- 


48.3 


0.35 


W+iets 


18.2 


2.7 ± 0.4 


ti 


16.4 


0.34 


WW 


28.0 


10.6 ± 0.1 


WZ 


19.2 


1.08 


Wy 


2.0 


0.03 


Predicted Backg 


round 1.19 x 10^ 


15.6 ± 0.4 


ZZ ^ e+n'+i 


2.9 


0.02 


ZZ e+Tuu 


8.9 


4.03 


Predicted Total 


1.19 X 10^ 


19.6 ± 0.4 


Data 


118,850 


28 


TABLE II: Number of predicted dimuon events and yield ob- 


served in data after the dilepton selection and after the re- 


quirement on The uncertainties in the final column are 


statistical only. 


If not present, the statistical uncertainty is 


negligible. 






Sample 


dilepton selection 


]^'rp requirement 


Z £+£- 


1.30 X lO'' 


0.1 ± 0.1 


Z T+T~ 


53.3 


0.09 


VK-I-Jets 




< 0.01 


tt 


16.0 


0.21 


WW 


32.0 


9.7 ± 0.1 


WZ 


18.3 


0.82 


Predicted Backg 


round 1.30 x lO'^ 


10.9 ± 0.3 


ZZ £+£-£'+£ 


2.89 


0.00 


ZZ ~* £+£-w 


9.48 


3.39 


Predicted Total 


1.30 X lO'^ 


14.3 ± 0.3 


Data 


127,960 


15 



The ZZ signal is separated from the remaining back- 
grounds with significant ^'-^ using a likelihood with the 
following input variables: the invariant mass of the dilep- 
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ton pair, Mu (for the dielectron channel), the prob- 
abihty resulting from a refit of the measured lepton mo- 
menta under the constraint that their dilepton mass gives 
the Z mass (for the dimuon channel), the transverse mo- 
mentum of the higher pT lepton, pr^^^ , the opening angle 
between the dilepton pair and the leading lepton, A0, 
and the cosine of the negative lepton scattering angle in 
the dilepton rest frame, cos(6'*). Figures [S] and [7] show 
the data and predicted distributions of the variables used 
in the likelihood for the dielectron and dimuon channels 
respectively. Fig. [5] shows the likelihood distributions for 
signal and backgrounds after all selection requirements. 
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FIG. 6: Distribution in the dielectron channel of the input 
variables of the likelihood discriminant for data and MC. In- 
variant mass of the dilepton system (a), pt of the leading 
lepton (b), the opening angle between the lead lepton and 
the dilepton system (c) , and the cosine of the scattering angle 
of the negative lepton in the dilepton rest frame (d) . 
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FIG. 7: Distribution in the dimuon channel of the input vari- 
ables of the likelihood discriminant for data and MC. prob- 
ability for the kinematic fit to the dilepton mass (a) , pT of the 
leading lepton (b), the opening angle between the lead lepton 
and the dilepton system(c), and the cosine of the scattering 
angle of the negative lepton in the dilepton rest frame (d). 
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FIG. 8: The likelihood distributions for signal and back- 
ground: dielectron events (a) and dimuon events (b). All 
selection requirements have been applied. 



SYSTEMATIC UNCERTAINTIES 



Systematic uncertainties are evaluated separately for 
the dielectron and the dimuon samples and for each of 
the data taking periods. The uncertainties affecting the 
overall scale factor of the MC cross sections are canceled 
out by normalizing to the data before the Ifi'j, cut. The 
remaining systematic uncertainties contributing to the 
significance and cross section are dominated by the nor- 
malization of the VF-l-jets background, the uncertainty 
on the WW cross section, the lepton resolution, and the 
number of Z events surviving the cut. The dominant 
uncertainties are listed in Table IIIII in which Az is the 
acceptance times efficiency for Z/7* — > and Azz 

is the acceptance times efficiency for ZZ/^* (.^l~vv 



where contributions from t^t^ decays are in- 

cluded. The large uncertainty on the VF-|-jets and re- 
maining Z are due to the uncertainties on the jet to lep- 
ton misidentification rate used in the normalization of 
the M^-|-jets background and the small statistics avail- 
able after the ^'j, cut (for both). Varying the parameters 
of the electron and muon smearing in the MC shows that 
the effect on the final result is within the statistical un- 
certainty in almost all bins. It is therefore propagated as 
an uncertainty in the shape of the likelihood, as are the 
contributions from jet energy resolution and the shape of 
the Z Z pt spectrum. 
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CROSS SECTION AND SIGNIFICANCE 

A negative log-likelihood ratio (LLR) test statistic is 
used to evaluate the significance of the result, taking as 
input the binned outputs of the dielectron and dimuon 
likelihood discriminants for each of the two data taking 
periods. A modified frequentist calculation is used [16| 
which returns the probability of the background only 
fluctuating to give the observed yield or higher (p-value) 
and the corresponding Gaussian equivalent significance. 
The combined dielectron and dimuon channels yield an 
observed significance of 2.6 standard deviations (2.0 ex- 
pected) , as reported in Table IIVI 

Because of the background normalization method de- 
scribed earlier, the measurement of the ZZ l^l^vv 
production cross section can be therefore expressed in 
terms of relative number of events with respect to the 
Z -> i+e.- sample. 

We define a background hypothesis to include the dis- 
tributions of the predicted backgrounds shown in TablesU 
and nil and a signal hypothesis to include these back- 
grounds and the events from the ZZ process. 

To determine the cross section the likelihood distri- 
bution in the data has been fitted allowing the signal 
normalization to float. The scale factor / with respect to 
the SM cross section and its uncertainty are determined 
by the fit. The ZZ production cross section is computed 
by scaling the number of events predicted by the MC to 
obtain that in data: 



a{ZZ) = a{Z) 



A, 



A 



ZZ 



Nz 



Azz is found to be 4.73 ± 0.03% in the dielectron chan- 
nel and 4.91 ± 0.03% in the dimuon channel. We as- 
sume the theoretical cross section for Zj^* — > t^t^ 
in the mass window 60 < M« < 130 GeV: (t(Z) = 
256.6^^2^0 Pb 0, [il- Using the ratio of the ZZ to 
ZZj-j* cross sections computed with MCFM [li| at NLO, 
we scale the ZZj^* cross section down by 3.4% to give 
a pure ZZ cross section. The resulting cross section for 
^ + A is 

a{ZZ) = 2.01 ± 0.93(siat.) ± 0.29(sys.) pb. 

This can be compared with the predicted SM cross sec- 
tion of 1.4 ± 0.1 pb [3 at = 1.96 TeV. 



TABLE III: The values assumed by the dominant systematic 
uncertainties for the various Monte Carlo signal and back- 
ground samples in the dielectron and dimuon channels. 



Systematic uncertainty 


dielectron 

(%) 


dimuon 

(%) 


-I- Jets normalization 


16 




vv vv ana vv zj 






Theoretical cross sections 


7 


7 


Number of Z events surviving 






the cut 


18 


3 


Systematic uncertainty on the 


Uncertainty 


cross section 


(%) 




theoretical cross section 


+2.0 
-5.0 


+2.0 
-5.0 


Az/Azz ratio from pdf 


uncertainties 


1.8 


1.8 


Az/Azz ratio from modeling 






of the veto efficiency 
Az/Azz ratio from modeling 


0.8 


0.8 


of the ZZ pt spectrum 


3.0 


3.0 



TABLE IV: Estimated significance for background only to 
fiuctuate to at least the observed yield for the combined di- 
electron and dimuon channels in the two data taking periods. 





expected (a) 


observed (a) 


p- value 


0.0244 


0.0042 


significance 


2.0 


2.6 



This is in agreement with the standard model prediction 
of 1.4 pb [11. 
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CONCLUSION 



We performed a measurement of the production cross 
section of ZZ t^i~vi? using 2.7 fb^^ of data col- 
lected by the DO experiment at a center of mass energy 
of 1.96 TeV. We observe a signal with a 2.6 standard de- 
viations significance (2.0 expected) and measure a cross 
section cr(pp ZZ) = 2.01 ± 0.93(stat.) ± 0.29(sys.) pb. 
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